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Abstract 

There has been an increasing interest on the concept of Inelastic Dark Matter 
(iDM) - motivated in part by some recent data. We describe the constraints 
on iDM from the results of the two phase dark matter detector ZEPLIN- 
II, which has demonstrated strong background discrimination capabilities 
(>98%). Using previously published estimates of the ZEPLIN-II residual 
background, the iDM limits presented here exclude a significantly larger iDM 
parameter space than the limits derived without background subtraction, 
suggesting the exclusion of iDM signal claims at >99% C.L., for Weakly 
Interacting Massive Particles (WIMPs) masses >100 GeV. 
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A liquid xenon (Xe) detector with two phase readout structure is well 
suited to the search for dark matter [H This is due in part to the high 
atomic mass of Xe and in part to the powerful background rejection (>98%) 
of the two phase detector, which was first applied to the dark matter search 
by the UCLA Torino group in 1996. We show here the key steps on how the 
background was measured for the ZEPLIN-II detector, as were first reported 
in Ref. j^, and applied the recent results to place constraints on Inelastic 
Dark Matter (iDM) model. 
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1. ZEPLIN-II Operation 



ZEPLIN-II is a dual pliase (gas and liquid) Xe detector designed specif- 
ically to identify nuclear recoils induced on Xe atoms by the Weakly Inter- 
acting Massive Particles (WIMPs); one of the popular candidates for dark 
matter in the universe [3], The detector is located in the UK Bou- 
bly underground laboratory which provides about 2800 m water-equalvent 
shielding for a 10^ reduction in cosmic muon flux [6|. The central detec- 
tor was constructed out of low-radioactive materials at UCLA, with under- 
ground ambient radioactivity carefully studied and modeled at Boulby 
Additional background rejection is afforded by the dual phase design which 
measures both scintillation and ionization components for each event in the 
liquid xenon. Since background electron recoils differ from the signal nuclear 
recoils in the ratio of these components by about a factor of 3, background 
rejection can be achieved at >98% efficiency ZEPLIN-II detector took 
data for 31.2 days with 7.2 kg fiducial mass and obtained the final data set 
reported in Ref. [3]. The WIMP signal acceptance box is defined between 
5 to 20 keV electron-equivalent (keVee), encompassing 50% of the total sig- 
nal region determined by prior calibration runs. 29 events were observed in 
the acceptance box at the end of first dark matter run. However, since the 
acceptance box is flanked on both sides by background events, it is to be 
expected that some of the events in the box come from the tail distributions 
of the background events. A careful analysis of the background distribu- 
tion was performed, and it was learned that 28.6 ± 4.3 background events 
were expected in the acceptance box. The detailed analysis and results are 
documented in Ref. [sj, and are summarized here. 



2. ZEPLIN-II Background Measurement 

The background events in the S2/S1~300 region above the acceptance 
box are 7 events, as learned from calibration runs with 7 source ^°Co and 
americium/beryllium (AmBe) neutron source. To quantify the event rate 
distribution, events with energy between 5 to 20 keVee were selected and 
binned according to its Log(S2/Sl)-k(E) values, where k(E) is the 50% nu- 
clear recoil acceptance value of log(S2/Sl) for the energy of each event. In 
another words, events with Log(S2/Sl)-k(E)<0 reside in the acceptance box. 
The 7 event rate distribution is well characterized by a Gaussian with an off- 
set, where the offset accounted for the coincidental events arising from the 
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high trigger rate (~70 Hz) encountered during cahbrations. These coinci- 
dental events were not expected to contribute to the background level in the 
acceptance box, since they were not observed in the science data, which oc- 
curred at a much lower trigger rate (~2 Hz). As such, only the Gaussian 
component of the fit is used in background estimation, where the Gaussian is 
scaled to the overall event count in the science data and integrated up to the 
50% nuclear recoil boundary. In total, 16.1 ± 2.9 7 events were estimated to 
be in the acceptance box between 5 to 20 keVee- 

The second background events found below the acceptance box were nuclear 
recoil events with diminished S2 originated from radon contamination in 
SAES getter, as confirmed by an dedicated Rn measurement. These recoils 
occurred near the charged PTFE walls which stripped some of the charges 
from the electron cloud that constituted S2. The diminished S2 subsequently 
lowered the position reconstruction accuracy, which resulted in a fraction of 
these events being wrongly placed within the fiducial volume. To quantify 
the number of misplaced events, a distribution of reconstructed radii of all 
events within the acceptance regions in (S2/Sl)-energy parameter space was 
studied. This distribution is well fitted with a Gaussian and is extrapolated 
into the fiducial volume. Finally the extrapolated distribution is integrated 
to give the expected number of radon-related events, which is 12.5 ± 2.3. 
Refer to Table [1] for a summary of background estimation. 



Energy Range 


Observed 


7-ray ^^Co) 


Rn-initiated 


Total 


5-10 keVee 


14 


4.2 ± 2.4 


10.2 ± 2.2 


14.4 ± 3.3 


10-20 keVee 


15 


11.9 ± 2.7 


2.3 ± 0.5 


14.2 ± 2.7 



Table 1: Overall expectation values in the nuclear recoil acceptance window compared to 
observed counts. The errors are derived directly from fit uncertainties. 



3. Inelastic Dark Matter Model 

The inelastic dark matter (iDM) [sl, [qI, i3] model explored the possibil 
ity of an altered kinematic of the WIMP-nucleus interaction in an attempt to 
reconcile the detection of annual modulation observation 
from all other experiments such as CDMS XENONIO 
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and null results 
and ZEPLIN- 



II jsj. iDM assumes two basic properties: a) WIMP x has an excited state 
X*, with a mass m^. — m-^ = 5 ~ 100 keV recoil-energy (keV^). b) Elastic 
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scatterings of x off target nucleus N are suppressed compared with the in 



elastic scattering xN— s>x*N. Refer to Ref. [lOj, [ij] and [l5| for discussions 
on particle models that may give rise to this scenario. 

Compared to the elastic case, the introduction of splitting 5 in the scatter- 
ing kinematics increases the minimum velocity (imin of X to scatter with a 
deposited recoil energy E/j. The new /^^m can be shown to be 



where m^r is the mass of the target nucleus and n is the reduced mass of the 
x/target nucleus system. For direct detection experiments, the key conse- 
quences of iDM are: 

1. Differential event rate is lowered since minimum velocity [3min is in- 
creased for a given recoil energy. 

2. Low-energy events are suppressed in the spectrum of events. In the 
elastic case, event rate is the highest in the low recoil energy, while the 
splitting 5 ~100 keV^ in iDM scenario can suppress or even eliminate 
low energy events. 

3. Annual modulation of signal is enhanced. In the extreme case, where 
the X particles only have enough minimum velocity in the summer but 
not in the winter, the modulation can be 100%, albeit at a significant 
lower overall event rate. 

4. Heavier target nuclei are even more favored. While it is true that in 
the elastic case heavier target nuclei gives rises to a higher differential 
event rate, inelastic scattering favors heavier target more since (3min is 
lowered significantly for heavier elements. 

These key characteristics are quantified in Fig [H An effective test to the 



WIMPs discovery signals reported in Ref. with Iodine (I) target (A=127) 
is to use a different target with comparable atomic mass. ZEPLIN-II with 
Xe (A=132) target is a good experiment for this test. 



4. ZEPLIN-II Inelastic Dark Matter Limits 

The first ZEPLIN-II dark matter run recorded 29 events in the accep- 
tance box with an expected background of 28.6 ± 4.3. The result is con- 
sistent with a null experiment with an upper limit of 10.4 using Feldman 
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Cousins method [16[. We maintained that the number of observed events 
in the acceptance box is consistent with the tail distributions from the two 
background populations. On the other hand, a confidence interval derived by 
taking all 29 events as iDM signals without background subtractions, such 



as the pMax method detailed in Ref. [17[, may not be the best approach to 



interpret ZEPLIN-II results since PMax is in general more suitable for ex- 
periments with small number of events where a statistical characterization 
of background is more challenging. We maintain that we understood our 
background accurately enough to use Feldman Cousins' method, which has 
an additional advantage of guaranteed correct coverage. 
Fig. [2] shows iDM ZEPLIN-II exclusion limits and DAMA signal contours 
for different WIMP masses. DAMA contours of 90% (red, Ax^=6.25) and 
jreen, Ax^=11.34) C.L. were constructed following the treatments by 
using iodine quenching factor (ql) of 0.085, where = Ege/ql. For 
ZEPLIN-II, the zero field quenching factor qXe is modified by the presence of 
the electric field, and the conversion equation becomes Eji = Eee/qXe-(Se/Sn), 
where the scintillation quenchirig of electron and nuclear recoils due to the 
electric field are Se=0.54±0.02 [3] and S„=0.95±0.05 Assuming Stan- 
dard Halo Model (SHM) parameters outlined in Ref. 




2ll and 22 , and a 



constant qXe of 0.19 (giving Er = Eee/0.36), the resultant ZEPLIN-II in- 
elastic dark matter limits are shown in blue solid lines, for the WIMP masses 
between 70 GeV < m^ < 250 GeV. Using these parameters, it can be seen 
that these limits constrained a larger iDM parameter space than those pre- 
viously reported in Ref. [sl, suggesting the exclusion of published claims for 
iDM signals at >99% C.L. for WIMP masses >100 GeV. Also shown in pink 
dotted lines are ZEPLIN-II limits assuming the the most recent energy de- 
pendent qXe values reported in Ref. [i^] , which estimated an increase in these 
values at high recoil energies (qXe ~0.24 at 50 keV^). The new qXe mea- 
surements lowered the exclusion power of ZEPLIN-II results, but at WIMPs 
masses >180 GeV there are still strong indications of exclusion of discov- 
ery claims at >99% C.L. Incidentally, uncertainty in ql, which is reported 
as ql=0.09±0.01 Q, has a significant effect on DAMA/LIBRA preferred 
regions in iDM case. In general, a higher ql would increase the discrepancy 
between DAMA/LIBRA and the other experiments, while a lower ql would 



lead to a better agreement. Refer to Ref. [IJ] for a more detailed discussion of 
the effects of ql, as well as other cosmological parameters, on the constraints 
of iDM. 



5 



5. Conclusion 

In conclusion, the published ZEPLIN-II background analysis reported in 
Ref. showed that the events observed in the acceptance box were statisti- 
cally consistent with the tail distributions from two background populations; 
namely the 7 and radon-induced events. As such it is natural to take the 
background subtraction approach to analyze the events in the acceptance re- 
gion, from which we declared ZEPLIN-Il a null experiment with a 90% CL. 
upper limit on WIMP signal number at 10.4. This translated to iDM exclu- 
sion limits that are significantly more stringent than the previously published 
values in Ref. j^, which adopted PMax method to derive an upper limit. Using 
parameters from Standard Halo Model, ZEPLIN-II limits suggest the exclu- 
sion of published claims for iDM signals at >99% C.L., for WIMP masses 
>100 GeV. 
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Xenon target, elastic and inelastic rates 
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Figure 1: (top) Energy spectrum of WIMP events on Xe target for both elastic scattering 
(red, dotted) and inelastic scattering (blue, line) models, (middle) Inelastic event rate 
suppression, compared to the elastic function of the splitting 5. (bottom) Xe 

annual modulation fraction as a function of S assuming Veart/i=227±14.4 km/s. The above 
plots assumed galactic escape velocity Vescape=500 km/s. 
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Figure 2: Preferred regions for DAMA/LIBRA results and ZEPLIN-II exclusion limits. 
The blue solid lines are ZEPLIN-II exclusion limits at 90% C.L., assuming xenon zero-field 
quenching factor qXe=0.19. The pink dotted lines are the exclusion limits assuming the 
most recent qXe reported in Ref. [20|. The DAMA/LIBRA contours are at 90%i[red) and 
99% (green) C.L. with standard ql=0.085, following the analysis done in Ref. [Sj. These 
plots assumed Wescape= 500 km/s. 
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